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Dipole Polarizability of 120 Sn and Nuclear Energy Density Functionals 
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The electric dipole strength distribution in 120 Sn between 5 and 22 MeV has been determined 
at RCNP Osaka from polarization transfer observables measured in proton inelastic scattering at 
Eo = 295 MeV and forward angles including 0°. Combined with photoabsorption data a highly 
precise electric dipole polarizability aD( 120 Sn) = 8.93(36) fm 3 is extracted. The dipole polarizability 
as isovector observable par excellence carries direct information on the nuclear symmetry energy 
and its density dependence. The correlation of the new value with the well established aD( 208 Pb) 
serves as a test of its prediction by nuclear energy density functionals (EDFs). Models based on 
modern Skyrme interactions describe the data fairly well while most calculations based on relativistic 
Hamiltonians cannot. 

PACS numbers: 21.10.Ky, 25.40.Ep, 21.60.Jz, 27.60.+j 


The nuclear equation of state (EOS) describing the en¬ 
ergy of nuclear matter as function of its density has wide 
impact on nuclear physics and astrophysics [l| as well 
as physics beyond the standard model [§, 3]. The EOS 
of symmetric nuclear matter with equal proton and neu¬ 
tron densities is well constrained from the ground state 
properties of finite nuclei, especially in the region of sat¬ 
uration density po ~ 0.16 fm -3 [4J]. However, the de¬ 
scription of astrophysical systems as, e.g., neutron stars 
requires knowledge of the EoS for asymmetric matter fi¬ 
ll] which is related to the leading isovector parameters 
of nuclear matter, viz. the symmetry energy (J) and its 
derivative with respect to density (L) [sj]. For a recent 
overview of experimental and theoretical studies of the 
symmetry energy see Ref. E3 ■ In spite of steady exten¬ 
sion of knowledge on exotic nuclei, just these isovector 
properties are poorly determined by fits to experimental 
ground state data because the valley of nuclear stability 
is still extremely narrow along isotopic chains EMI. 


Thus one needs observables in finite nuclei specifically 
sensitive to isovector properties to better confine J and 
L. There are two such observables, the neutron skin r s ki n 
in nuclei with large neutron excess and the (static) dipole 
polarizability qd. 


The neutron skin thickness r s ki n = {r) n — ( r) v defined 
as the difference of the neutron and proton root-mean- 
square radii (r) n ,p is determined by the interplay between 
the surface tension and the pressure of excess neutrons on 


El. Studies within nuclear 
show for all EDFs a strong 
correlation betwe en r s ki n and the isovector symmetry en- 

The most studied case so far is 


the core described by L 
density-funtional theory 
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ergy parameters 
208 Pb, where r s ki n has been derived from coherent pho¬ 
toproduction of 7T° mesons |2(| , antiproton annihilation 


2ll . |23], proton elastic scattering at 650 MeV |23| and 
295 MeV [3], and from the dipole polarizability [23 • A 
nearly model-independent determination of the neutron 
skin is possible by measuring the weak form factor of nu¬ 
clei with parity-violating elastic electron scattering |26| . 
Such an experiment has been performed for 208 Pb but 
the statistical uncertainties are still too large for serious 
constraints of the neutron skin [27]. 

A particularly useful experimental observable to con¬ 
strain the large theoretical uncertainties on J and L is 
ao [H which can determined by a weighted integral over 
the photoabsorption cross section er a b s @ 

Tic fcr^ hs 8 tt /'dH(El) 

“ D = 275 / = y / -g—M 
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where E x is the excitation energy and B( El) the reduced 
electric dipole transition strength. It is the aim of this 
letter to present a new experimental result for ao in a 
heavy nucleus, 120 Sn. This data point is then used to¬ 
gether with the well established CKD( 2 CKS Pb) to scrutinize 
EDFs. 

The El response is dominated by excitation of the 
isovector giant dipole resonance (IVGDR) well known in 
many nuclei from photoabsorption experiments. Because 
of the inverse energy weighting in Eq. ©, «d also de¬ 
pends on the low-energy strength studied mainly with the 
( 7 , 7 ') reaction. However, extraction of the El strength 
from ( 7 , 7 ') data is rather model-dependent SiOi . 

Recently, polarized inelastic proton scattering at 295 
MeV and at forward angles including 0° has been es¬ 
tablished as a new method to extract the complete El 
strength in heavy nuclei from low excitation energy across 
the giant resonance region 25|. In this particular kine¬ 
matics selective excitation of El and spin-Ml dipole 
modes is observed. Their contributions to the cross sec¬ 
tions can be separated either by a multipole decomposi¬ 
tion analysis (MDA) j^J or independently by measure¬ 
ment of a combination of polarization transfer observ¬ 
ables (PTA) (2fj. Good agreement of both methods was 
demonstrated for the reference case 208 Pb where values of 
r s kin and L derived from aD( 208 Pb) conform with results 
from other methods [HJ. 

All EDFs agree on showing strong correlations between 
a d, t skin? J, and L, but the actual predictions of qd for 
given J and L values differ considerably. While the result 
for 208 Pb [H] already excluded many older Skyrme in¬ 
teractions, modern Skyrme-Hartree-Fock (SHF) and rel¬ 
ativistic mean-field (RMF) models can be brought into 
agreement e.g. by changing J, which can be varied over 
a certain range without deteriorating the fit of the inter¬ 
action parameters 0. Experimental information on qd 
in other nuclei is therefore essential to further constrain 
the isovector part of the EDF interaction. We note that 
some information on El strength distributions in heavy 
neutron-rich nuclei is available but an extraction 

of qd from these results is completely model-dependent, 
in contrast to the data discussed here. 

Here, we report on a measurement of the electric 
dipole response in 120 Sn with polarized proton scatter¬ 
ing based on a PTA covering excitation energies 5 — 22 
MeV. El strength in 120 Sn below 5 MeV was measured 
by [2LH ') 132] and above neutron threshold by ( 7 , xn) 
|38|-|40( experiments. A combination of all available data 
enables a precise determination of an- The El strength 
has also been determined from a MDA of the (p,p') cross 


sections [4l| but photoabsorption cross sections had to be 
included as constraints and therefore the result - in con¬ 
trast to the PTA - is not independent from these data. 

The experiment was performed at the RING cyclotron 
facility of the Research Center for Nuclear Physics 
(RCNP), Osaka University, Japan. Details of the exper¬ 



FIG. 1. (Color online), (a) Double differential cross sections 
(black) and decomposition into non-spinflip (red) and spinflip 
(blue) parts of the 120 Sn(p, p r ) reaction at Eq = 295 MeV and 
9 — 0° — 2.5°. The green line shows the cross sections due to 
excitation of the ISGQR estimated as described in the text, 
(b) Total spin transfer from Eq. @. 


imental technique can be found in Ref. [42j] . A polarized 
proton beam was accelerated to Eq = 295 MeV and scat¬ 
tered protons were momentum-analyzed with the Grand 
Raiden spectrometer [43| placed at 0° covering an angu¬ 
lar and excitation energy range of 0° — 2.5° and 5 — 22 
MeV, respectively. A 11 isotopically enriched (98.4%) self- 
supporting 120 Sn foil with a thickness of 6.5 mg/cm 2 
served as a target. The beam intensity was 1 — 2 nA 
with an average polarization of 0.7. 

A decomposition of spinflip and non-spinflip cross sec¬ 
tions can be achieved |44| by the combined information 
of the polarization transfer observables Dll, Dss and 
Dnn |45( determined in a secondary scattering experi¬ 
ment. Since Dss an d Dnn are indistinguishable at 0°, 
only Dll and Dss were measured in the present experi¬ 
ment. It is convenient to introduce the total spin transfer 


E = 


3 — 2 Dss — Dll 


( 2 ) 


which takes values of zero for non-spinflip and one for 
spinflip transitions. Because of the different reaction 
mechanism these can be identified with El (Coulomb 
excitation) or Ml (spin-isospinflip part of the proton- 
nucleus interaction) excitations, respectively. 

Figure [ID a) displays the measured cross sections (black 
circles) in 400 keV bins. The bump structure centered at 
E* ~ 15 MeV corresponds to the IVGDR. The extracted 
total spin transfer [Fig. [T](b)] is almost zero in this en¬ 
ergy region as expected for Coulomb excitation and ap- 
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proaches maximum values of about 0.2 around 9 MeV 
(the location of the Ml spinflip resonance [46J) and above 
18 MeV. The decomposition into non-spinflip and spinflip 
parts is shown in FigQJa) by red and blue circles, respec¬ 
tively. The non-spin-flip cross sections contain an small 
E2 contribution (green line) from nuclear excitation of 
the isoscalar giant quadrupole resonance (ISGQR). It was 
determined using the isoscalar B( E2) strength distribu¬ 
tion [Zzj as described in Ref. [4lJ and never exceeds 4% 
in a single bin. 



FIG. 2. (Color online). Comparison of the -B(El) strength 
distribution in 120 Sn determined by the present work (red 
circles) and in ( 7 , 2 : 71 ) experiments (blue circles [38]], green 
circles [39}], and black circles 0 )- 


Figure [2] shows the R(El) strength distribution (red 
circles) deduced from the AS = 0 cross sections assuming 
semiclassical Coulomb excitation [48]. The photoabsorp¬ 
tion data converted to B(E1) strength are shown as blue 
[38 j. green j®. and black [hj circles, respectively. All 
data agree well with each other. Near the IVGDR maxi¬ 
mum slightly smaller values are found in the present work 
but they still accord within the experimental uncertain¬ 
ties. Since the MDA analysis below neutron threshold 
does not depend on the photoabsorption data, the 1?(E1) 
strength up to 9 MeV can be compared to the present 
result. Good agreement of the summed strengths from 
MDA [48.7(29) e 2 fm 2 ] and PTA [54.1(41) e 2 fm 2 ] and the 


corresponding contribution to an is observed. 

We now turn to the determination of the electric dipole 
polarizability in 120 Sn taking into account all available 
data. The energy region below 5 MeV makes a negligible 
(< 0.1%) contribution to an [37|. Results for the energy 
region from 5 to 10 MeV are taken from the present work 
|49} and amount to 1.12(7) fm 3 , contributing about 12.5% 


to the total value. The main contribution, 7.00(29) fm 3 , 
stems from the IVGDR region, where the present results 
and those from Refs. 38 40( were averaged between 10 
and 22 MeV. Between 22 and 28.9 MeV data are available 
from Ref. 0], 0.51(6) fm 3 . Finally, the polarizability at 
even higher energies up to 135 MeV was taken from a 
nat Sn( 7 ,xn) experiment [50] neglecting an isotopic de¬ 
pendence. The contribution, 0.31(10) fm 3 , is small but 
non-negligible considering the final precision achieved. In 
total, we find aD( 120 Sn) = 8.93(36) fm 3 , where the er¬ 
ror contains the statistical and systematic uncertainties 
of all data used. 

Having now at hand precise data for OD( 120 Sn) and 
Q ,D( 208 pb), we use them to scrutinize the performance 
of a broad variety of EDFs from SHF and RMF, all val¬ 
ues including pairing at the BCS level. The theoretical 
an values are computed from the static response to an 
external dipole field. Figure [3] displays the EDF results 
for aD( 208 Pb) versus OD( 120 Sn) together with the exper¬ 
imental values indicated by yellow bands. 

Panel (a) collects SHF results for a couple of widely 
used parametrizations (see caption). Although taken 
from very different sources, all SHF results together show 
a strong correlation between the theoretical q;d values in¬ 
dicated by the dashed black line. The actual position on 
the line is determined by J and the large span of results 
along the line demonstrates the uncertainty in J. Note, 
however, that the majority of SHF results resides nicely 
within the experimentally allowed yellow square. The 
outliers are all rather old parametrizations adjusted be¬ 
fore appearance of the many data on neutron rich nuclei. 
The fact that the linear trend goes right through the ex¬ 
perimentally allowed square and that most parametriza- 
tion lies within indicates that the isovector density de¬ 
pendence of SHF is realistic. 

Two parametrizations (SV-min, RD-min) are shown 
together with error bars from statistical analysis El®- 
These are larger than the experimental uncertainties 
demonstrating that the data provide indeed useful con¬ 
straints on the isovector parameters [l 2]. 

Simple error bars hide the linear correlation dis¬ 
cussed above. This can be better visualized by se¬ 
ries of par ametrizations with systematically varied J 
0, 13, 49]. We do this in connection with RMF ap¬ 
proaches shown in panel (b). Unlike SHF, there is greater 
variance in modeling density dependence for RMF. We 
consider three variants thereof: the density dependent 
point-coupling model (DD-PC) ,61], the density depen¬ 
dent meson-exchange model (DD-ME) [62], and non¬ 
linear meson coupling in FSU [63[. For all three cases we 
show series with varied J (open symbols) and the best 
fit (full symbols with error bars where available). For 
better comparability, the series DD-ME, DD-PC, and 
RD were fitted to the same data pool as SV-min and 
RD-min [0. ( Fit p rocedures for the FSU family are de¬ 
scribed in Refs. [63j,|64[). However, the actual fit strategy 
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FIG. 3. (Color online). Correlation of the experimental a® 
values for 120 Sn and 208 Pb with uncertainties shown as yel¬ 
low bands, (a) Comparison with models based on Skyrme 
interactions SkM* [5l| (red square), SkP [52j] (blue triangle), 
SkT6 [Hij (red diamond), SG-II [541 (blue circle), SkI3 [5a] 
(green triangle), SLy6 [56j] (red triangle), BSk4 [5711 (green 
circle), SV-bas (blue diamond) and UNEDF2 [58l] (blue 
square). The SV-min [ll| (red circle) and RD-min [59| 
(green square) interactions additionally provide theoretical 
error bars [3, [6til | . The dashed black line indicates the cor¬ 
relation between both qd values, (b) Comparison with rel¬ 
ativistic mean held models DD-PC-min (blue squares) and 
DD-ME-min (black circles), both from Ref. [l2|], FSU [63j] 
(red diamonds), and FSU2 [64[ (green triangle). Full symbols 
denote the results of optimum parameter sets. Open symbols 
show results varying the symmetry energy parameter J. The 
dashed lines serve to guide the eye. The dashed black line 
from (a) is repeated for direct comparison of SHF and RMF 
models. 


seems to be of lower importance as the original standard 
parametrizations DD-ME2 [62| and PC-1 6l| lie again on 
the corresponding lines in the plot. All sets are strongly 
correlated with nearly linear trends, however, with dif¬ 
ferent offset depending on the form of the EDF. While 
the SHF series goes approximately through the center 
of the experimental square, all RMF chains are off cen¬ 
ter, two of them just touching the square. Only DD-PC 
comes closer and only one DD-PC parametrization (with 
J = 32 MeV) lies within the correlation box. The best- 


fit parametrizations (full symbols) are all outside. This 
indicates that RMF models still need to be improved in 
the isovector channel [l2j , although the modern, density- 
dependent functionals already constitute large progress 
in this respect in comparison to older RMF functionals 



J_I_l 

8 9 10 

a D ( 120 Sn) (fni) 


FIG. 4. (Color online). Relationship between a® and r s ki n 
for 120 Sn predicted by the SV-min [l3j (red circle) and RD- 
min (59j (green square) interactions. Full symbols are the 
results of the optimum parameter sets and open symbols cor¬ 
respond to a variation of the symmetry energy parameter J 
as in Fig. [3[b). Dashed lines are to guide the eye. The hori¬ 
zontal lines denote the range of r s ki n values compatible with 
the experimental polarizability shown as yellow band. 


Using the strong correlation between a® and r s kin (HI 
one can derive the neutron skin thickness of 120 Sn from 
EDFs capable to describe the data in Fig. [3] A similar 
analysis has been performed for 208 Pb [18(. Since the 
models are not independent, rather than averaging (as 
done in Ref. [01) we take the SV-min (red circles) and 
RD-min (green squares) results as representative and es¬ 
timate the theoretical uncertainties. Figure [4] shows the 
predictions of the correlation between r s kin and crp- As 
for the relativistic models, a variation of J (open sym¬ 
bols) is compatible with the optimum fits (full symbols). 
The range of values consistent with the experimental po¬ 
larizability indicated by the horizontal lines corresponds 
to r s kin = 0.148(34) fm. The result is in good agree¬ 
ment with values extracted from measurements of the 
spin-dipole resonance (66j, 0.18(7) fm, and proton elastic 
scattering [67| , 0.16(3) fm, while antiproton annihilation 
[68} finds a much smaller value, 0.08(+3)(-4) fm. 

In summary, we have measured polarized proton in¬ 
elastic scattering off 120 Sn at very forward angles and 
extracted the El strength distribution between 5 and 
22 MeV by an analysis of polarization transfer observ- 
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ables. Combining the present results with ( 7 , xn) data, 
the dipole polarizability could be extracted with a pre¬ 
cision of 4%. The correlation with the polarizability of 
208 Pb [25| provides an important test of EDFs indispens¬ 
able for the extraction of properties of the symmetry en¬ 
ergy in neutron-rich matter. While modern Skyrme inter¬ 
actions can describe the data, in contrast to most RMF 
calculations. With the typical theoretical uncertainties 
indicated, the combined data from 208 Pb and 120 Sn pro¬ 
vide an important constraint to improve the description 
of static isovector properties in EDFs. 


Considering the importance of polarizability data, a 
systematic study at different shell closures and explo¬ 
ration of the role of deformation is called for. One im¬ 
portant future project is a systematic measurement of «d 
covering the range of stable tin isotopes Together 

with a new measurement of relativistic Coulomb excita¬ 
tion of the neutron-rich tin isotopes 124-134 Sn at GSI [70] 
a unique set of data will be available to investigate the 
impact of neutron excess on the formation of a neutron 
skin in a set of nuclei with similar underlying structure. 
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